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1. Introduction 

A flavor changing neutral current (FCNC) process is an apparent transition between quarks 
of different flavor but equal charge. In the standard model (SM), these processes occur through 
higher order loop diagrams. The SM branching fraction for such FCNC decays is quite small. It is 
interesting to study such rare decays because the contribution from physics beyond the SM may be 
sizable. In order to be able to observe rare heavy flavor decays it is essential to produce a sufficient 
number of B hadrons. At the Tevatron pp collider, various types of B hadrons can be produced 
with large production cross sections of ^(100 jib) [jl|]. However, the inelastic pp cross section is 
10 3 times larger than that of B hadron production. Interesting events have to be extracted from a 
high track multiplicity environment. Therefore, detectors need to have very good tracking, vertex 
resolution, wide acceptance, good particle identification and highly selective triggers. 

In this paper, we report on recent results from studies on rare B decays via FCNC processes, 
B°,s /i + jU~, B + -)• K+n+jj,-, B° K*°n + jj,- and fi° ->• 0jU+jU", performed by the CDF and 
DO collaborations. 

2. B° {s) ^ii+ii- 

The FCNC B decays with dimuon in the final state are highly suppressed by the helicity factor, 
(m^/mfi) 2 . In the SM, the branching fraction of 5° ->■ ji + ll~ decays is (3.6 ±0.3) x 10~ 9 [@]. 
B° — > decay is further suppressed by the ratio of CKM matrix elements, |\4</Atj| 2 @ |J, 

giving a branching fraction in the SM of (1.1 ±0.1) x 10~ 10 [|0. The decay amplitude can be 
enhanced by several orders of magnitude in some extensions of the SM. Since the predicted rate 
for these processes in the SM is beyond the current experimental sensitivity at the Tevatron, the 
observation of these decays would necessarily imply physics beyond the SM. Improved limits on 
the branching fraction can be used to set limits on the parameter space of supersymmetric models 
and other new theories. It is also important to measure the ratio of the branching fractions between 
5° — > jli + ;U~ and B° — > jU + /i~ to test the Minimal Flavor Violation [0]. 

CDF has performed a search for 2?9\ — > jJ. + H~ with 3.7 fb _1 of data []|]. Their acceptance 
is increased compared to the previous analysis [[J, in addition to increased luminosity, by using 
events where muon candidates cross the mid-plane of the central outer tracker (COT) where plastic 
inserts to maintain wire spacing create a zone with lower trigger efficiency. This additional trigger 
acceptance increases the total trigger acceptance by 12%. Backgrounds from hadrons misidentified 
as muons are now suppressed by selecting muon candidates using a likelihood function. This func- 
tion tests the consistency of electromagnetic and hadronic energy with that expected for minimum 
ionizing particles and the differences between extrapolated track trajectories and muon system hits. 
In addition, backgrounds from kaons that penetrate through the calorimeter to the muon system, 
or decay in flight outside the COT, are further suppressed by a loose selection based on the mea- 
surement of the ionization per unit path length, dE / dx. After applying their baseline selections, a 
Neural Network (NN) multivariate classifier is constructed using the following variables to achieve 
further separation of signal from background: (1) the measured proper decay time, (2) the proper 
decay time divided by the estimated uncertainty, (3) the 3D opening angle between the dimuon 
flight direction vector and the displacement vector between the interaction point and the dimuon 
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vertex, (4) the fi-candidate track isolation, (5) the B-candidate transverse momentum, (6) and the 
transverse momentum of the lower momentum muon candidate. The NN is trained using back- 
ground events sampled from the sideband regions and signal events generated with a CDF simula- 
tion. The NN output distributions of 5° signal and sideband background events are shown in Fig. []]. 
The search is performed in a two dimensional grid in dimuon mass and NN space to improve 
the sensitivity. The expected background is obtained by summing contributions from the combina- 
torial continuum and from B — > h + h'~ decays (h + and h'~ represent a charged kaon or pion), which 
peak in the signal invariant mass region and do not occur in the sidebands. The contribution from 
other heavy-flavor decays is negligible. The combinatorial background is estimated by linearly ex- 
trapolating from the sideband region to the signal region. The B — > h + h'~ contributions are about a 
factor of ten smaller than the combinatorial background. The dimuon invariant mass distributions 
for three different NN output ranges are shown in Fig. [|. The limit on the branching fraction is 
computed by normalizing to the number of reconstructed B + — > J/yK + events. A sample of about 
19,700 B + — > J/\i/K + events is collected to serve as a normalization channel using the same base- 
line requirements. CDF extracts new upper limits 3$(B° S — > jU + jU~) < 4.3 x 10~ 8 (3.6 x 10~ 8 ) and 
&{B Q -> jU + AO < 7.6 x 10~ 9 (6.0 x 10~ 9 ) at the 95% (90%) C.L., which are currently the world's 
best upper limits for both processes. 

DO has also performed a search for B® — > using 6.1 fb _1 of data The event selec- 

tion criteria are similar to those of CDF. The muon selection has been updated with respect to the 
previous analysis yielding 10% higher acceptance while keeping the fraction of misidentified 
muons below 0.5%. To further suppress the background, the following discriminating variables are 
used: (a) the transverse momentum of the B ( s ' candidate, (b) the 3D opening angle, (c) the transverse 
decay length significance, (d) the decay vertex fit % 2 , (e) the smaller impact parameter significance 
of the two muons, (f) and the smaller transverse momentum of the two muons. A Bayesian Neural 



Network (BNN) [H, |lfj multivariate classifier with the above variables is constructed to distinguish 



signal events from background. The BNN is trained using background events sampled from the 
sideband regions and simulated signal events. The distributions of the BNN output j8 for the B ( s ' 
signal and the sideband events as well as the B + — > J j\\fK + control sample are shown in Fig. ||. 
The 5° -> \l + \L- signal region is defined to be 0.9 < j3 < 1.0 and 5.0 GeV < m m < 5.8 GeV. 
Two-dimensional (2D) histograms of mna vs. j8, dividing the signal region into several bins, are 
prepared to improve the sensitivity relative to using a single bin. The dominant source of back- 
ground dimuon events is from decays of heavy flavor hadrons in bb or cc production. To study 
this background contribution, inclusive dimuon Monte Carlo samples with PYTHIA [ |TT| ] generic 
QCD processes are generated. The dimuon background events can be categorized by two types: 

(i) B(D) — > n + vX,B(D) — > pL~vX' double semileptonic decays where the two muons originate 
from different b(c) quarks, yielding dimuon masses distributed over the entire signal region, and 

(ii) B — > n + vD,D — > n~vX sequential semileptonic decays, resulting in m w predominantly be- 
low the B hadron mass. The simulated dimuon mass distributions for both background sources are 
parametrized using an exponential function to estimate the number of background events in the sig- 
nal region after fitting the dimuon mass in the data sideband regions in each /3 bin. The uncertainty 
on this background estimate is dominated by the statistical uncertainty of the sideband sample 
(10% ~ 35%). The B — >■ h + h'~ background contribution is negligible. The limit on the branch- 
ing fraction is computed by normalizing to the number of reconstructed B + — > J /yK + events. A 
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Figure 1: CDF data. Left: Distributions of the NN output for simulated B { ] — > M + M signal and sideband 
events in data. Right: The dimuon invariant mass distribution for events satisfying all selection criteria for 
the final three NN output ranges. 



sample of about 46,800 B + — > J j\\fK + events is selected using the same baseline selections. The 
observed distributions of dimuon events in the highest sensitivity region are shown in Fig. [2[ The 
observed number of events is consistent with the background expectations. The resulting limit is 
@[B° S -)• jU + M") < 5.1 x 10~ 8 (4.2 x 10~ 8 ) at the 95% (90%) C.L. The simulated mass resolution 
of the DO detector for the B® — > /i + jU~ is ss 120 MeV and is therefore insufficient to readily sepa- 
rate B*' from B° leptonic decays. In this analysis, it is assumed that there are no contributions from 
B° — > jU + /i~ decays, since this decay is suppressed by |V f rf/V f i| 2 ~ 0.04. 

Aside from the background uncertainty, the largest uncertainty of 15% common to CDF and 
DO B*' — > A t+ jU _ analyses comes from the fragmentation ratio between B + and B® [12]. Both 
CDF and DO observe no evidence of physics beyond the SM in B ^ — > /i + jU~ decays, but they 
significantly improve the existing limits providing tighter constraints for the parameter space of 
several SM extensions. 



3. 5° -r K*°ji + ji-, B+ -> K+n + n and 5° <^/i + M 

Another sensitive probe of non-SM physics accessible at the Tevatron is to measure the FCNC 
b —> stl semileptonic decays. Although there is no helicity factor in these decays, the SM branch- 
ing fraction is still small, of order ^"(lO" 6 ), and there could a sizable contribution to the decay 
rate from physics beyond the SM. In addition, the kinematic distributions could be modified. The 
lepton forward-backward asymmetry and the differential branching fraction as a function of dilep- 
ton invariant mass in the decays B — > K*U differ from the SM expectations in various extended 
models. The amplitudes of the effective Wilson coefficients Cj, C9, and C10 may interfere with 
the contributions from non-SM particles, where C7, C9, and Ciq represent the amplitudes from 
the electromagnetic penguin, the vector electroweak, and the axial-vector electroweak contribu- 
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Figure 2: DO data. Left: Distributions of j3 : (a) B s — » signal and sideband events, (b)B + -^J/\\/K + 

data and simulation. Right: The distribution of in the highest sensitivity j5 region (a), and the distribu- 
tion of j3 in the highest sensitivity region (b) for data (dots with uncertainties), expected background 
distribution (solid line), and the SM signal distribution multiplied by a factor of 100 (dotted-dashed line). 
The dimuon background contributions from the B(D) — > jj, + vX ,B(D) — > jj,~vX' decays (dashed line) and 
the B -> n + vD,D -)• n~vX decays (dotted line) are also shown. 



tions, respectively. Among many b — > sjx + n~ decays, the exclusive channels B + K + n + n and 
B° — > K*°}i + ii~ have been observed and studied at Belle and BaBar At the Tevatron it 
is also possible to search for analogous decays — > 0jU + jU~ and A° — > AjU + jU~. 

CDF reports new results on measurements of B + — > K + n + n~, B° — > K* Q pi + }i~ and B° s — > 
tyjl^ ll~ decays JT5|]. Event selection starts from selecting two oppositely charged muon candi- 
dates with a transverse momentum greater than 1.5 GeV or 2.0 GeV depending on the trigger 
selection. B — > hpL + pL~ candidates are then reconstructed by adding h = K + , K*°, or (j) candi- 
dates to the dimuon pair. The K*° is reconstructed in the mode K*° — > K + n~ and the is re- 
constructed as — >■ K + K . The kaon and pion candidates are required to be consistent with the 
time-of-flight (TOF) and dE jdx combined log-likelihood probability of each particle hypothesis. 
The muon candidates are further purified by the muon-likelihood. To enhance separation of sig- 
nal from background, an Artificial Neural Network multivariate classifier is constructed. Figure || 
shows the invariant mass distribution for each rare decay. By performing an unbinned maximum 
log-likelihood fit of the B invariant mass distribution, 120± 16 events for B + — > K + pL + pL , 101 ± 12 
events for B° — > K* ii + ji~ , and 27 ± 6 events for B® — > are found, with 8.5a, 9.7a and 

6.3a statistical significance, respectively. This is the first observation of tP s — > 0/x + /i _ decay. Us- 
ing the corresponding B — > J j\\fh modes as a reference, CDF determines the following absolute 
branching fractions: SB(B + K + \i + \x-) = [0.38 ±0.05 ±0.03] x 10~ 6 , m{B° ->■ ^*V + AO = 
[1.06 ±0.14 ±0.09] x 10~ 6 and ^(5° -> 0/x + /i _ ) = [1.44 ± 0.33 ±0.46] x 10~ 6 , where the first 
uncertainty is statistical and the second is systematic. The differential branching fractions with 
respect to q 2 = ml^c 2 for B° — > K* Q ji + }i~ and B + — > K + ji + n~ are also calculated as shown in 
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^ CDF Run II Preliminary L=4.4fb 
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Figure 3: CDF data. The B invai-iant mass of B° -> K*°fi + ^,- (left), B + -> A"+ju + jU" (middle) and 
5° (right). 
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Figure 4: CDF data. Differential branching fraction of B° K*°fi + ^~ (left) and B + -> K + pL + ii~ (right). 
The hatched regions are vetoed to prevent contributions from charmonium. The solid lines are the SM 
expectation, which use maximum- and minimum- allowed form factor. The dashed line is the averaged 
theoretical curve in each q 2 bin. 



Fig. |j. The forward-backward asymmetry (Apg) and K* Q longitudinal polarization (F£) are ex- 
tracted from cos 6^ and cos 6k distributions, respectively, where 6^ is the helicity angle between 
/i + (jU — ) direction and the direction opposite to the B (B) direction in the dimuon rest-frame, and 
6k is the angle between the kaon direction and the direction opposite to the B meson in the K*° rest 
frame. The fit results of Fl and Aps for B° — > K* pL + pL~ are shown in Fig. ||, as well as Apb for 
B + — > K + /j. + jj.~. Both Fl and Apb are consistent with the SM expectation, and also with an exam- 
ple of a SUSY model [|T^]. These results are also consistent and competitive with the B-factories' 
measurements rtT3l[T3]. 
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Figure 5: CDF data. F L (left) and A FB (middle) fit results as a function of q 2 for B° -> K*°n + n~, and 
A/rg (right) as a function of g 2 for B + — > K + jj. + fi~. The points show data. Solid (dotted) curve is the SM 
(an example of SUSY) expectation jig]. The dashed line is the averaged expectation in each q 2 bin. The 
hatched regions are vetoed to prevent contributions from charmonium. 



4. Conclusion 

CDF reports the first observation of the decay 5° — > 0jU + /i~, results on the forward-backward 
asymmetry measurements in fi° — > K* Q }1 + }1~ and B + — ► K + ii + }i~ decays using 4.4 fb _1 of inte- 
grated luminosity, and upper limits on the branching fraction £$(B° S — > jj, + jj,~)< 4.3 x 10~ 8 (3.6 x 
1(T 8 ), m{B° -»• jU+jU") < 7.6 x 1(T 9 (6.0 x at the 95% (90%) C.L. using 3.7 ftT 1 of inte- 

grated luminosity. DO also reports a new upper limit ,^{B° S — y /i + jU _ )< 5.1 x 10 _8 (4.2 x 10~ 8 ) 
at the 95% (90%) C.L. using 6.1 fb _1 of integrated luminosity. There is not a significant discrep- 
ancy from the SM expectation in the FCNC B decays so far. All these results are done with about 
half or less the statistics expected in a year from now and the analyses are continuously improved, 
therefore we could expect large improvements in the near future. 

References 

[1] CDF Collaboration, D. Acosta et al, Phys. Rev. D 71, 032001 (2005) [arXiv:hep-ex/0412071]. 

[2] A. J. Buras, Prog. Theor. Phys. 122, 145 (2009) [arXiv:hep-ph/0904.4917]. 

[3] N. Cabibbo, Phys. Rev. Lett. 10, 531 (1963). 

[4] M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 49, 652 (1973). 

[5] CDF Collaboration, Note 9892. 

[6] CDF Collaboration, T. Aaltonen et al, Phys. Rev. Lett. 100, 101802 (2008) 
[arXiv:hep-ex/0712.1708]. 

[7] DO Collaboration, V. M. Abazov et al, Phys. Lett. B 693, 539 (2010) [arXiv: 1006.3469 [hep-ex]]. 

[8] DO Collaboration, V. M. Abazov et al, Phys. Rev. D 76, 092001 (2007) [arXiv:hep-ex/0707.3997]. 

[9] R. M. Neal, Bayesian Learning of Neural Networks (Springer- Verlag, New York, 1996). 

[10] P. C. Bhat and H. B. Prosper, "Bayesian Neural Networks" in Statistical Problems in Particle Physics, 
Astrophysics and Cosmology, ed. L. Lyons and M. K. Unel, (Imperial College Press, London, 2006). 



7 



Rare B decays at the Tevatron 



Masato Aoki 



[11] T. Sjostrand et al, Comput. Phys. Commun. 135, 238 (2001). 
[12] W.-M. Yao et al, Journal of Physics G 33, 1 (2006). 

[13] BELLE Collaboration, J. T. Wei et al, Phys. Rev. Lett. 103, 171801 (2009) [arXiv:hep-ex/0904.0770]. 
[14] BABAR Collaboration, B. Aubert et al, Phys. Rev. D 79, 031 102 (2009) [arXiv:hep-ex/0804.4412]. 
[15] CDF Collaboration, Note 10047. 

[16] A. Ali, P. Ball, L. T. Handoko and G. Hiller, Phys. Rev. D 61, 074024 (2000) [arXiv:hep-ph/9910221]. 



8 



